Role of saddles in mean-field dynamics above the glass transition 
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fSJ . Recent numerical developments in the study of glassy systems have shown that it is possible to 

» I ' give a purely geometric interpretation of the dynamic glass transition by considering the properties 

CIh, of unstable saddle points of the energy. Here we further develop this program in the context of a 

mean-field model, by analytically studying the properties of the closest saddle point to an equilibrium 
configuration of the system. We prove that when the glass transition is approached the energy of 
the closest saddle goes to the threshold energy, defined as the energy level below which the degree of 
instability of the typical stationary points vanishes. Moreover, we show that the distance between a 
typical equilibrium configuration and the closest saddle is always very small and that, surprisingly, 

■ it is almost independent of the temperature. 

o 

The glass transition occurs when the relaxation time of a substance increases upon cooling of many orders of 
magnitude in a very narrow interval of temperature, without the onset of any crystalline order. Even though dramatic 
changes in the mechanical properties of the sample occur, it is impossible to define a strict transition temperature, 
because the dynamic process leading to the glassy phase is continuous, albeit very sudden. This is indeed one of the 
^ ■ most tricky points in the study of glassy systems: the so-called glass transition cannot actually be associated to the 
genuine divergence of any dynamic or thermodynamic quantity. In fact, its very definition as a reference temperature 
is based upon convention: in supercooled liquids, it has been agreed to fix the glass transition temperature Tg at the 
point where the viscosity of the sample (that is one of the macroscopic manifestation of the relaxation time) is of 
^ ; order lO^^P. 

■ On the other hand, in the case of fragile super-cooled liquids there is at least another important value of the 
O ' temperature which is useful to describe and interpret experimental data, that is the temperature Tmct where Mode 

, . Coupling Theory (MCT) locates a purely dynamic transition |^ . Such a transition is spurious, since what is observed 
in real experiments and simulations is just a dynamical crossover from a diffusive regime to an Arrhenius (or super 
y—i Arrhenius) one. However, MCT describes well the dynamics of fragile liquids for T > Tmct and, even though the 
^ ' MCT transition is smeared out in reality, still the MCT temperature remains a meaningful reference value marking 
t I the border between purely diffusive and activated slow dynamics 

■ The lack of a strict dynamic transition is not common to all glassy systems. It has been discovered in the past 
^ ' Q that some mean-field models for spin-glasses display a phenomenology quite similar to the one of real structural 

I glasses and supercooled liquids, but for a notable difference: in these mean- field systems there is a true divergence 
^-H ' of the relaxation time, with no associated thermodynamic anomaly. This fact makes the definition of a dynamic 
, critical temperature Td completely unambiguous for these models. Besides, the dynamical equations which describe 
' the behaviour of these systems above coincide with those obtained by MCT Thus, for these models MCT is 
^ I exact, and Td therefore coincides with Tmct- The most deeply studied among these systems is the p-spin spherical 
^ ' model l^-[l^, henceforth indicated as pSM. 

A key feature of the pSM is the possibility to explain the dynamic glass transition at Td as the result of a purely 
geometric transition taking place in the energy landscape of the system |^-[To|]. At a given energy density, called 
threshold energy Etu, there is a qualitative change in the stability properties of the landscape: below Eth and down to 
Q , the ground state energy Eq minima dominate, whereas above Eth unstable saddles are the most numerous stationary 
points of the Hamiltonian. It can be proved that in such a system the dynamic glass transition occurs when the 
equilibrium energy density becomes equal to the energy density of the threshold states at that temperature In 
other words, in the pSM the dynamic glass transition at Td and the geometric transition at Eth are essentially two 
\ faces of the same phenomenon. Due to this fact, the structure and properties of unstable stationary points in the 
pSM have been the object of a number of investigations in recent years | p^ , p^ , p^ . 

Infinite lifetime metastable states cannot exist in finite dimensional systems, and therefore we cannot expect to find 
a divergence of the relaxation time with no associated thermodynamic transition in non-mean-field models. However, a 
strict geometric transition at a threshold energy may very well occur also in more realistic systems, such as supercooled 
liquids, even if its dynamical counterpart is smeared out by the finite dimensional nature of the system (i.e. by the 
finite lifetime of the threshold minima). If this were true, the dynamic crossover at Tmct, which in supercooled 
liquids marks the onset of activated glassy dynamics, would actually be the manifestation of a more fundamental and 
sharply defined geometric transition occurring at a certain critical threshold energy. 
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This scenario has been numerically investi gate d very recently inj^,|l6| for various Lennard- Jones (LJ) systems, 
and from a more speculative point of view in [|l7| . In particular, in |15f| a well defined threshold potential energy has 
been located and associated to the onset of glassy dynamics in the system. Moreover, it has been directly shown in 
1^ that the energy landscape of LJ models and that of the pSM are indeed very similar. These studies seem therefore 
to confirm the idea that, even in realistic systems, the dynamic crossover observed at Tmct is the consequence of the 
sharp change in the topology of the energy landscape at the threshold energy. 

The approaches of JT^I-pTt all have as a vital starting point the assumption that the time evolution of the system in 
the phase space is in some way infiuenced by the nearby saddle points of the potential energy. However, this fact has 
not been directly proved, and the circumstantial evidences are mainly of a numerical nature. In particular, it may 
be objected that the trajectory of the system at equilibrium is never even close to saddles, especially above the glass 
transition, where it may be argued that free diffusion in the phase space implies that the energy landscape and its 
stationary points are completely irrelevant (for the relevance of saddles in zero temperature dynamics see [p^). Even 
in the pSM there was up to now little evidence of any direct connection between dynamics of the system equilibrated 
above Td and saddle points of the energy above Eth (see, however, the approach of Furthermore, even assuming 

that the dynamic trajectory stays somewhat close to saddles, it remains to be directly demonstrated that such objects 
do play a role in the transition. More precisely, one should prove that the properties of these supposedly close saddles 
indeed display some anomaly at the dynamic transition. 

The aim of the present work is therefore to analytically investigate what is the role of saddles in the equilibrium 
dynamics of the pSM above T^- In order to do this we will introduce a tool, which allows for the exact location of 
the closest saddle points to an equilibrium configuration at temperature T. In this way we will be able to study how 
the properties of these closest saddles vary with the temperature when the dynamic transition is approached, thus 
answering some of the questions raised above. 

The Hamiltonian of the pSM is given by, 

N , N 



H 



J^.■■■^,n,■■■T,^+0{\/N) , (1) 
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where the spins satisfy the spherical constraint — ^- The quenched couplings Ji^ - i^ are Gaussian distributed 

random variables with variance = p\/2N'P^^ . By means of the Lagrange method, we can find the stationary points 
of the Hamiltonian on the sphere and therefore write the equations satisfied by the saddle points of H with energy 
density i?, 

1 ^ 

- Y ■Jk,^,■■■^,n, ■■■n^^ETk^O , fc = 1, . . . , TV . (2) 

^' i2---ip 

In general, the number Af{E) of solutions of equations (|^) with energy density E is exponentially large in the size of 
the system N. Thus, the quantity which is normally computed is the complexity (or configurational entropy), defined 
as the logarithmic density of this number, S(i?) — log JV{E). The nature of the saddle points of H is in principle 
not only specified by their energy density E, but also by their instability index K, that is the number of negative 
eigenvalues of the Hessian matrix. However, previous studies of the pSM have shown that there is a well defined 
relation between energy density and index K{E), and that at any fixed enerm level E only stationary points with 
index K{E) dominates the energy landscape in the thermodynamic limit P, [l^ , p^ . Therefore, by fixing the energy 
density of a saddle point to E, we are automatically fixing its index to K{E). A crucial feature of the pSM is that 
the typical saddles index is extensive, K = 0{N), as long as the energy density is above a value called threshold, Eth, 
while K = for E < Eth- This means that minima dominate over saddles below the threshold, while saddles of index 
K{E) > are the most numerous stationary points for E > Eth- In this sense, we can say that at Eth a geometric 
transition takes place. More precisely, if we introduce the index density k = K/N , we have. 
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Note that k{E) is a monotonically increasing function of the energy E. Remarkably, when the equilibrium energy 
density of the system becomes equal to the internal energy density of the threshold minima the system undergoes a 
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dynamic glass transition, which we will indicate with T^. To better specify this statement, we have to distinguish 
between the bare energy ii^ of a minimum, and its internal energy U{T), that is the energy of a system equilibrated 
in that minimum at temperature T Q Of course, the quantity U(T) is equal to the bare energy E plus a vibrational 
contribution due to thermal fluctuations. At the dynamic glass transition Td we have that Ueq{Td) — UthiTd), where 
Ueq(T) is the global equilibrium energy density of the system, and Uth{Td) — Eth + vibrati 



ons. 



Our aim in this paper is to analyze the structure of the saddle points around an equilibrium configuration thcrmalizcd 
at temperature T > Td- We therefore need a notion of distance to give a meaning to this statement. Given two 
configurations a and r we define a co-distance, or overlap, qa-r, as, 

1 ^ 

i 

Similar configurations have g ~ 1, while different ones have g ~ 0. Our strategy will be to fix a reference equilibrium 
configuration a and compute the complexity of the saddles points r close to it as a function of their overlap qar- 
The value of the overlap where this quantity goes to zero will give the distance of the closest stationary points to a. 
Indeed, for larger overlaps, i.e. smaller distances, a negative complexity indicates a vanishing probability of finding a 
stationary point. 

In order to do this we have to calculate how many saddles r, with a given energy E, happen to have an overlap q 
with a reference equilibrium configuration a. Clearly, this number formally depends on a itself and on the disorder J . 
However, as always done in similar calculations po| , pll[l]}| , we can assume that in the thermodynamic limit N ^ oo 
this quantity is self-averaging with respect to the distribution of a and J, and therefore we can average it over the 
Gaussian distribution of the disorder (indicated with a bar) and over the equilibrium distribution of a at temperature 
T. In this way we can define the saddles complexity as. 



S,(g,i?,/3)^ 1 J 
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(4) 



with, 



Z{P) = dae 



-l3H(a) 



and where integration is carried out over spherical configurations only. To understand equation (M) it is convenient 
to read it right-to-left: first, under the r integral, we calculate using the standard method of [g2| the number of 
solutions of equations (j^), putting an extra constraint on the overlap they must have with cr. Second, we take the 
logarithm of this quantity, and we average it over the equilibrium distribution of a. Finally, we average everything 
over the disorder J. As we can see, depends on the temperature T = 1//? at which the reference configuration a 
is equilibrated, on the energy E of the saddles r we are counting, and finally on the overlap q between a and r. We 
stress that, by construction, a is an independent equilibrium configuration, irrespective of the energy and the distance 
of the saddle t. 

In order to perform the averages in it is convenient to use the replica method, writing 

Z-^ = lim 

(log(.)) = lim llog((.n . (5) 

rn^O m 



In this way we have. 



^ In other terms, if rrii is the local magnetization of the system equilibrated in the minimum, and q = 1/N^^m1 is the 
self-overlap (i.e. the magnetization norm, which is related to thermal fluctuations), we can define a bare magnetization as 
rhi — The bare energy density is then given by _E = -^H{rhi). 
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with a = 1, . . . , n and a = 1, . . . , m. The exphcit calculation of from equation (^) can be performed by using the 
standard tools of the replica method: variational parameters are introduced and the integrals are evaluated exactly 
in the limit ^ 00 by means of the steepest descent method. Of course, this is possible thanks to the mean-field 
nature of the model. Here, we will skip most of the details and just state the final result. The interested reader may 
refer to |jri|, where a technically similar calculation is performed. The full expression for the saddle complexity is: 

S,(q,£;,/3) = 1 + + i log 
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with, 



fli = {xi - xo){l - r) + {yi - y2f 

^2 = (a;o + -r) + [yi ~ yo)[2(2/o - qw) - [r - q^){yi - yo)/il - r)]. 



(8) 

As customary in the context of the replica method, the set of variational parameters x = [xq, xi,yo,yi, w, r) is fixed 
be means of the steepest descent equations 9I]s/i9x = 0, which we have solved numerically. We remark that the 
expression above for Es is only valid in the regime E > Eth- In showing the results we will assume p = 3. 

First of all, we are interested in studying the behaviour of T,s{q, E, T) as a function of g, at fixed E and T. In this 
way we can define an overlap qi^{E, T) where Eg goes to zero: this overlap gives the distance of the closest saddle with 
energy E to an equilibrium configuration at temperature T. In Figure 1 we plot Eg as a function of the overlap q, for 
T — Tfi and E = Et^ . At this temperature many properties of equilibrium landscape are known and an interpretation 
of the results is therefore much simpler. At Td the system equilibrates inside a threshold state with bare energy 
density Eth and self-overlap (i.e. largeness) qth (|]- In naive terms we can then imagine that our typical equilibrium 
configuration a lies in a well whose largeness is given by qth and whose bottom is at energy density Eth- In this case 
it is evident that the closest stationary point to a is precisely the bottom of the well. The point where the complexity 
goes to zero must therefore give the overlap between the center of the threshold minimum r and one of its typical 
equilibrium configurations a. This overlap can be easily computed by noticing that q^r — '^ii'^iTi) — J4 'T^i'Tii 
where the thermal average is restricted to a threshold state, and indicates the local magnetization of that state. 
In the pSM this local magnetization can be expressed directly in terms of the well minimum as — ~^qthTi 1^ (see 
also footnote 1) and we immediately get g^r = J qth- Consistently with this result we find. 



qQ{Eth,Td) 



(9) 



This result can be appreciated in Figure 1. 
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FIG. 1. Saddles complexity as a function of the overlap q, for E — Eth and T — Td- Inset: same curve in linear-log scale. 
The complexity goes to zero at go ~ x/Qth ~ 0.71, for p = 3. 
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Moreover, a careful analysis of S^. for q ^ qo shows that, 



^M,Eth,Td)^iq~qo)' , q-^qo. (10) 

Note that the exponent is the same as found in for the approach to zero of the constrained complexity of threshold 
Thouless- Anderson-Palmer (TAP) solutions. This is a consistency check for the present calculation. 

In the light of our original aim to find the closest saddles, it is interesting to plot the value of go a-s a function of 
the energy density E of the saddles we are counting. We expect this curve to have a maximum at a value Es(T) 
corresponding to the energy of the closest, and thus the most relevant, saddles. Accordingly, in Figure 2 we plot 
qa{E) as a function of E for different values of the temperature. In particular, the full curve represents qo{E) for 
T = Tci: as we see, this is a steadily decreasing curve having its maximum at the threshold energy. This means that, 
as previously said, the closest stationary point to an equilibrium configuration at the glass transition is a threshold 
minimum, Es{Td) = Eth- However, if we now increase the temperature T of the equilibrium configuration, we expect 
the energy of the closest saddle to increase as well, together with its instability index. In other words, the higher the 
temperature of a, the higher will be the energy, and thus the degree of instability, of the closest possible saddle (we 
remind that K{E) is a monotonic increasing function of E, with K{Eth) — 0). This hypothesis is confirmed in Figure 
2: the maximum of these curves moves to the right as the temperature is increased, disclosing a well defined relation 
Es(T), which we will analyze carefully later. For the moment, let us note that a further consistency check of our 
calculation is that for T = cx3 we find Es — (see Figure 2). At very high temperatures the equilibrium configuration 
a is just a random configuration of the system, therefore the closest stationary points to it will be the most numerous 
ones in absolute terms. In the pSM it can be proved that the most numerous saddles have E = and K = N/2 
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FIG. 2. The overlap go where the saddles complexity goes to zero as a function of the energy E of the saddles, at five different 
values of the temperature above Td- The full vertical line marks the value of the threshold energy Eth- 



It is interesting to study the behaviour of the overlap qs where the curves qo (E) have their maximum. This overlap 
is a measure of the closeness of an equilibrium configuration a to its nearest saddle point r. What is surprising of 
Figure 2, is that by varying the temperature the value of q^ is almost constant. To better investigate this point we 
plot in Figure 3 gs as a function of p. We can see that qs is practically always constant, but for /3 ~ Pa, where it 
sharply jumps to y/qth- This fact means that the distance between an equilibrium configuration and its closest saddles 
is almost independent of the temperature. This value of the overlap is qs ~ 0.68, which is indeed quite high, being 
comparable to the overlap between equilibrium configurations and bottom of the minima below the glass transition 
(see, for example. Figure 1). This result answers one of the main questions raised in the introduction: above the 
dynamic glass transition, the equilibrium trajectory indeed stays always very close to unstable stationary points of 
the Hamiltonian, exactly as below Td it stays close to stable minima. On the other hand, it is clear from Figure 2 
that for T > Td the closest minima (i.e. the ones with E = Eth) are very far from the dynamic trajectory. In this 
sense, it is justified to say that the equilibrium dynamics of the system above the glass transition may be described 
as an evolution among the neighborhoods of saddle points, rather than among basins of the minima |17|| . 
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FIG. 3. The overlap qa of the closest saddles to an equilibrium configuration at temperature T, as a function oi P — 1/T. 
Inset: enlargement of the same plot for (i ^ (3d- 



The distance of the closest saddle does not change with T, but as we have seen, the energy density Eg does. In 
Figure 4 we plot this energy as a function of the inverse temperature and compare it with the equilibrium energy 
density of the system above the glass transition, i.e. Ueq{T) = —13/2 First of all we note that the energy of the 
saddle is always smaller than the energy of the equilibrium configuration, despite the two objects being so close in 
the phase space. This fact has been already noted in the context of a numerical study of a Lennard-Jones system in 
| p6[ . It is tempting to interpret the difference Ueq{T) — EsiT) as a pseudo- vibrational contribution of saddles, due 
to the fact that, even though X > 0, the largest part of the Hessian eigenvalues is positive, as long as Eth < E < 0. 
Unstable saddles are not trapping object, of course, but they may have a substantially long life-time provided that K 
is small enough. This phenomenon is at the basis of the pseudo-vibrational contribution of saddles. 

This last hypothesis is supported by another interesting result we find, that is, 

EsiT) ^ Eth , T^Td. (11) 

Therefore, the energy density of the closest saddles to an equilibrium configuration goes to the threshold energy 
density at the dynamic glass transition. Clearly, at the difference between Es{Td) and Ueq{Td) is given by the 
vibrational contribution of thermal fluctuations inside threshold minima, which is of order kgT. When T > 'we see 
that the two curves continuously approach one another, as the saddles instability index K increases with the energy. 
These results seem thus to suggest that, even though the system is not confined into any given saddle point, the 
disproportion between trapping and un-trapping directions, that is the fact that K < N/2, is sufficient to produce a 
vibrational contribution that we may broadly interpret as thermal fluctuations around saddles point. 




0.5 1 1.5 R 2 2.5 

FIG. 4. The bare energy Es of the closest saddles compared to the equilibrium energy Ueq- Inset: enlargement of the same 
plot for P ^ Pd- The slope of the curve changes in this regime. 
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A further support of this idea comes from the comparison of our results with the approach developed in , where a 
dynamic description of the pSM based on the concept of quasi-states is introduced. More specifically, for temperatures 
slightly above the dynamic transition Td, quasi-states are related to critical points of the TAP free energy (for a more 
detailed definition see jl^) and it can be shown that their global contribution gives rise to the paramagnetic free 
energy of the system. Interestingly enough, we found that the bare energy density of these quasi-states (defined as in 
footnote 1) is to a good degree of accuracy (within 1%) equal to the energy density Eg of the closest saddles. This 
indicates that the quasi-states introduced in p4| may be interpreted as our closest saddles plus the pseudo- vibrational 
contribution mentioned before. This interpretation confirms the idea, outlined in fl^ and made more explicit in the 
present paper, that above Td the paramagnetic state is made up of disjoint quasi-states around saddles. Equilibrium 
dynamics can thus be thought as evolution from one quasi-state (i.e. one saddle and its own neighborhood) to another 
one. A difference between the present approach and the one of ||l^ is that the pseudo-states of can be defined 
only very close to Td, while, as we have seen, closest saddles exist at any temperature, although, of course, we do not 
expect them to have any relevance for T ^ Td- 

Given the relation k{E) between index and energy of the typical saddles, we can introduce a temperature-dependent 
index by using the energy density Es{T) of the closest saddles, namely k{T) ~ k{Es{T)). Clearly, this index vanishes 
at the dynamic transition temperature T^, which is just another way of describing the geometric transition occurring 
at the threshold energy Eth- Close to Td we find that Eg is linear in (3 and this, together with the analytic form of 
k{E), implies that, 

k{T) cx (T - Tdf" , T^Td. (12) 

Summarizing, in the context of the pSM we have calculated the complexity of the saddle points at fixed overlap with 
a reference equilibrium configuration, above the dynamic glass transition. In this way we were able to identify what 
are the energy and distance of the closest saddles at any given temperature. We found that the distance between 
equilibrium configuration and closest saddle is almost independent of the temperature and is very small. Moreover, the 
energy of the closest saddles intersects the threshold energy at the dynamic glass transition. Finally, we interpreted 
the difference between equilibrium energy and energy of the closest saddle as a pseudo-vibrational contribution due 
to the fact that a number of trapping directions larger than the number of non-trapping ones may give rise to thermal 
fiuctuations around unstable saddle points. The present study supports the idea that dynamics in glassy systems for 
T > Td can be described in terms of evolution in the phase space among the neighborhoods of unstable saddles and 
strengthens the hypothesis that the glass transition, even in finite dimensional systems, is just the manifestation of 
the topological transition between saddles and minima dominated regions of the phase space. 
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[1] C.A. Angell, 'Strong and fragile liquids', in Relaxation in Complex Systems, K. Ngai and G.B. Wright eds., National 

Technical Information Series, Springfield VA, (1985). 
[2] U. Bengtzelius, W. G:otze and A. Sj:olander, J. Chem. Phys. 17, 5915 (1984); E. Leutheusser, Phy. Rev. A 29, 2765 

(1984). 

[3] C.A. Angell, J. Phys. Chem. Sohds 49, 1924 (1988). 

[4] T.R. Kirkpatrick and D. Thirumalai, Phys. Rev. B 36, 5388 (1987); T.R. Kirkpatrick and P.G. Wolynes, Phys. Rev. B 36, 
8552 (1987); Phys. Rev. A 35, 3072 (1987). 

[5] J. P. Bouchaud, L.F. Cughandolo, J. Kurchan and M. Mezard, Physica A 226, 246 (1996). 

[6] A. Crisanti and H.-J. Sommers, Z. Phys. B 87, 341 (1992). 

[7] A. Crisanti, H. Horner and H.-J. Sommers, Z. Phys. B 92, 257 (1993). 

[8] L.F. Cughandolo and J. Kurchan, Phys. Rev. Lett. 71, 173 (1993). 

[9] J. Kurchan, G. Parisi and M.A. Virasoro, J. Phys. I France 3, 1819 (1993). 
[10] A. Crisanti and H.-J. Sommers, J. Phys. I France 5, 805 (1995). 
[11] A. Cavagna, I. Giardina, and G. Parisi, J. Phys. A: Math. Gen. 30, 7021 (1997). 
[12] A. Cavagna, I. Giardina and G. Parisi, Phys. Rev. B 57 11251, (1998). 
[13] J. Kurchan and L. Laloux, J. Phys. A: Math. Gen. 29, 1929 (1996). 
[14] S. Franz and M.A. Virasoro, J. Phys. A: Math, and Gen. 33, 891 (2000). 

[15] K. Broderix, K. K. Bhattacharya, A. Cavagna, A. Zippelius, and I. Giardina, Phys. Rev. Lett. 85, 5360 (2000). 



7 



[16] L. Angelani, R. Di Leonardo, G. Ruocco, A. Scala, and F. Sciortino, Phys. Rev. Lett. 85, 5356 (2000). 

[17] A. Cavagna, Europhys. Lett. 53, 490 (2001). 

[18] A. Cavagna, J.R Garrahan L Giardina, Phys. Rev. B 61 3960, (2000). 

[19] K. Broderix, K. K. Bhattacharya, A. Cavagna, A. Zippelius, and L Giaxdina, in Disordered and Complex Systems, AIP 

Conference Proceedings 553, 23 (2001). 
[20] S. Franz and G. Parisi, J. Phys. I France 5, 1401 (1995). 
[21] A. Cavagna, L Giardina, and G. Parisi, J. Phys. A 30, 4449 (1997). 
[22] A.J. Bray and M.A. Moore, J. Phys. C 13, L469 (1980). 



8 



